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We present an experimental and theoretical study of electron tunnelling through quantum dots 
which focuses the attention on the amplitude of the current resonances as a function of magnetic 
field. We demonstrate that the amplitudes of the resonances in the tunnelling spectra show a 
dramatically different behaviour as a function of the magnetic field, depending on the angular 
momentum of the dot state through which tunnelling occurs. This investigation allows us to probe 
the details of the confined wave functions of the quantum dot. 



Quantum dots (QDs) are structures where electrons 
are confined in all directions. Ever since it has become 
possible to manufacture such structures there has been 
a considerable effort to develop spectroscopic tools capa- 
ble of detecting the energy spectrum of QDs. Among the 
different techniques used to probe the structure of zero- 
dimensional systems, magneto-transport is one of the 
most common because the presence of a magnetic field 
introduces a structure-independent quantisation that is 
superimposed on the one produced by the fabrication 
process 0], In all these studies the emphasis has been 
placed on the shift of the I-V features produced by the 
QD states as a function of the magnetic field. This type of 
analysis has been performed to probe the QD spectrum, 
but it is insensitive to the details of the QD wave func- 
tions. One notable exception has been the study by P.H. 
Beton et al. where the one dimensional wave functions 
were probed by means of the magneto-tunnelling. 

In this letter we will perform an experimental and the- 
oretical analysis of magneto-transport in vertical double- 
barrier systems enclosing a QD. We shall focus our at- 
tention on the amplitude of the current resonances as- 
sociated with the QD levels and on their variation with 
the applied magnetic field. We have observed a strik- 
ing difference in the behaviour of the amplitude of the 
current features. This is associated with the different 
angular momentum of the QD states through which the 
tunnelling occurs. 

In order to interpret this effect, we have developed a 
theoretical model which enables us to calculate the I- 
V characteristics of QD structures in the presence of a 
magnetic field analytically (^]. The model also shows 
that, in the presence of the magnetic field, a selection rule 
applies which severely limits or suppresses altogether the 
current flowing through the QD states with an angular 
momentum co- linear to B. 

The structures considered in this work consist of 
two 3D contacts with a QD inserted in the middle of 



a GaAlAs-GaAs-GaAlAs double barrier resonant tun- 
nelling heterostructure ( see figures 0a and |b). The 
fabrication of the samples is a combination of ion- 
implantation through metallic masks (used to define the 
quantum boxes lateral dimensions) patterned by electron 
beam and of epitaxial regrowth of the 3D GaAs top elec- 
trode 0. The GaAs quantum well has a width of 5.1nm; 
the well is embedded between two j4^o.33Gao.67^4s barri- 
ers of width 8.7nm. Si-doped contact layers with Si to 
2 x 10 17 cm -3 are formed on either side of the barriers. 
Care has been taken to prevent impurities inclusion in 
the quantum well region during the growth: the doped 
regions are separated from the barriers by a 20nm un- 
doped spacer layer. These layers have been grown in 
two steps: the first lOnm at a temperature of 540 °C to 
prevent Si segregation and the following lOnm plus the 
double barrier structure at 630 °C. 

In the following we will describe the results obtained 
on a nominal 35nm-radius dot at a temperature of 35mA" 
and with a magnetic field applied along the vertical di- 
rection of the device, i.e. along the direction in which the 
current flows. The main resonance current peak which is 
observed in larger devices is replaced by current plateaus 
(»,B,A), whose voltage positions and current amplitudes 
vary with B (see fig. |j) 

Before we further discuss these plateaus, we briefly fo- 
cus on the anomalously low threshold voltage (13mV), 
which is well below the expected value of 200mV. In 
fig. 0c we show the effect of LED illumination on the volt- 
age position of the first current resonance at T = 35mA. 
The TV characteristics have been recorded after each 
LED pulse. We clearly observe a persistent shift of the 
whole spectrum to lower voltages, which can be removed 
by increasing the temperature j^j. We believe that this 
behaviour is related to the ion-implantation process used 
to define the QDs which induces the formation of deep 
charged defects under the lateral barriers. As in con- 
ventional 2DEG heterojunctions, the light changes the 
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charge state of the deep defects, hence inducing a shift 
of the pinning of the Fermi level, lowering the threshold 
voltage. Despite the presence of such deep defects in the 
lateral barrier, we will show that the data presented here 
are unambiguously related to lateral confinement of the 
QD states. 

Let us first focus on the magnetic field dependence of 
the voltage positions of the current plateaus. In fig. ||a we 
show the dl/dV curves corresponding to the I-V spectra: 
the differential conductance traces reveal another reso- 
nance at high bias (T) in addition to the three main 
conductance peaks. Minor peaks labeled with a 'V are 
temperature independent and thus they are not related 
to QD states but to Local Density of State Fluctuations 
in the contacts [[§. 

The magnetic field displacement of the four conduc- 
tance peaks with B is compared in fig. ||b with the single 
particle spectrum of the QD where the energy eigenvalues 
are given by the Darwin-Fock formula: 

E nd jn d = huj d (2n d + \m d \ + 1) + \hjj c m d + E (1) 



where w d — y lu^ + ^ , lu is the strength of the size 
confinement, oj c is the cyclotron frequency and Eq is the 
energy of the quantum well ground state in the absence 
of lateral confinement. 

The best fit for the position of the peaks on the volt- 
age scale as a function of B gives Hlvq = IbmeV , with the 
voltage-to-energy conversion coefficient a = 0.5. The ra- 
dius of the first state is equal to r ~ ^j2h/m*ujQ ~ 12nm, 
which is in agreement with the nominal size of the QD 
taking into account the real lateral confinement potential 
profile induced by the ion-implantation inter-mixing pro- 
cess under the mask (tJ . Although single donors in GaAs 
have a comparable radius (lOnm), their excited states 
do not behave like the data of fig. |j|b. Coulomb blockade 
is also irrelevant because the distance between the two 
first peaks decreases with B instead of slightly increasing 
since B confines more the electrons, hence increasing the 
Coulomb repulsion. Thus we are dealing with features 
releated to QD states for which a parabolic confinement 
and a single-electron picture are good approximations, 
given the strength of the size confinement in our struc- 
tures ||. Spin effects are also negleted here since we do 
not observe such effects in our experiments. 

Let us focus now on the magnetic field dependence 
of the plateaus' current amplitudes. The fig. || shows 
that the amplitude of the first one decreases as B is in- 
creased, whereas that of the second one increases and 
then decreases. The third resonance amplitude decreases 
regularly. The current amplitudes of each QD state are 
calculated more precisely by integrating the conductance 
peaks in fig. [|a. The resulting curves have been reported 
on fig. |4^l for the first three main conductance peaks. In 
order to explain their behaviour, we have developed a 



method to calculate the I-V characteristics as a function 
of B [||. The QD eigenfunctions are: 

e im d tp 

$ dot (r, ip, z) = —j=-R nd>md (r)^ d {z) (2) 

where R nd .m d {f) is a function that can be expressed in 
terms of the hypergeometric function and *& d (z) is 

the ground state wave function of the one dimensional 
quantum well in the vertical direction. Similarly the ap- 
plied magnetic field produces highly degenerate ID Lan- 
dau subbands in the contacts. The spectrum of the con- 
tacts is given by: 

h 2 k 2 

En e ,m e ,k x = fuv e (2n e + m e + \m e \ + 1) + - — J (3) 

2m 

where lu £ = lu c /2 and k z is the momentum in the growth 
direction while the eigenfunctions are plane waves along 
the vertical direction and for the lateral part they are of 
the same shape as those of the QD. The only difference 
is that in the contacts only the magnetic field produces a 
quantisation whilst in the dot the effect of the magnetic 
field is superimposed on that of the size confinement. No- 
tice that the lowest Landau subband of energy E = huj e 
does not present strictly positive values of m e , the sec- 
ond subband of lateral energy E — Zhuj e contains only 
the value m e < 1, and in general the n e -th subband con- 
tains a maximum positive m = n e . 

Whenever the energy of a QD state falls between the 
Fermi level and the conduction band bottom, tunnelling 
through such a state becomes energetically possible and 
current starts flowing through it. The amplitude of the 
current that flows is determined by the overlap between 
the QD wavefunctions and the electron wavefunctions in 
the contacts which are at the same energy. The states in 
the contacts are highly degenerate ID Landau subbands 
which arc classified by the principal quantum number n e 
and the angular momentum m e along the direction of 
the field. The QD states are also classified in terms of a 
principal quantum number n d and the angular momen- 
tum m d . It follows that the tunnelling process must obey 
a selection rule, i.e. tunnelling is only possible between 
states with the same angular momentum J3j. Therefore, 
of the many degenerate subbands contained in each Lan- 
dau level in the contacts, only the one with the right value 
of angular momentum will contribute to the tunnelling 
through each QD state. 

The overlap integrals (Tl]] between the QD and the con- 
tacts' states can be evaluated analytically and for the 
lateral directions are given by: 

Ml d , m „ e = ^r( 7 )A-»-»-T(A-fc d r 

x (A - k e r<F(-n d , -n e , 7> ^ _ ^ _ - } )8„ d (4) 
where a a - 1 f (l m l+"^) ! \\ „, - I I 

wneie a d ,e - i+Im] i2\ m \n de l\m\l 2 ' > Ud ' e ~ V 2m*w ed (B)> 
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F is the hypergeometric function j|, A = 

ij e = „ ^ and 7 = \m\ + 1 (hereafter m = to^ and 

n = rid for clarity). This formula can be combined 
with well known expressions for the overlap integral in 
the vertical direction J|,[l2| to obtain the total overlap 
between each dot state and the contacts and therefore 
calculate the current within a transfer matrix approach. 
This method is well known to correctly reproduce the 
physics even in the resonant tunnelling regime of double- 
barrier systems. 

This technique has been adopted to calculate the I- 
V characteristics of 3D-0D-3D structures in the presence 
of an external magnetic field applied in the vertical di- 
rection, i.e. parallel to the current. Fig. ||b shows the 
behaviour of the current amplitude for different states. 
States n = and m = —2, — 1, 0, 1 are represented, cal- 
culated with the values of the parameters in accordance 
with those extracted from the fit using equation |l| of the 
experimental data. The amplitude for all states shows os- 
cillations as a function of magnetic field. Some of them 
are marked by an open triangle. These oscillations have 
the characteristic -g dependence and are due to the os- 
cillations in the density of states at the Fermi energy in 
the contacts as the magnetic field is changed. 

As the magnetic field is increased, the resonances show 
a rather different behaviour in terms of their amplitude. 
The resonance associated with the ground state shows 
an amplitude that decreases systematically as a func- 
tion of the magnetic field, while the first excited state 
n = 0,r7i = —1 initially increases in amplitude, at- 
tains a maximum, and eventually decreases. The state 
n = 0, m = —2 has the same qualitative features as the 
11 = 0,m = —1 one, but the maximum is more pro- 
nounced and shifted at higher magnetic field. 

The presence of the magnetic field produces two effects. 
Firstly it reduces the momentum along the direction of 
the current for the electrons in the contacts by moving 
the subbands' minima to higher energies. Secondly it 
confines the electrons towards the centre of the dot. The 
first effect tends to diminish the current, while the sec- 
ond increases it by increasing the overlap in the lateral 
direction between the QD wavefunctions and those of the 
contacts. Thus each Landau subband in the emitter gives 
a current which attains a maximum at a finite magnetic 
field. The position of this maximum depends on the rel- 
ative strength of the two effects. Moreover increasing the 
angular momentum m decreases the overlap integral and 
shifts this maximum to higher field. In the inset of fig. [|b 
we show the contribution of the first Landau level in the 
emitter for m = 0, —1, —2 to clearly illustrate this shift of 
the maximum. The total current is obtained by summing 
the contribution of all the occupied Landau levels in the 
emitter and it follows that the shift of the maximum cur- 
rent value from the origin B=0 T is a direct measurement 
of the magnitude of the angular momentum m. 



Perhaps even more striking is the totally different be- 
haviour of the states n = 0,m = 1 and n = 0,777 = — I 
when B increases. At B = these states are degener- 
ate and the radial part of the wave function is identical 
for both of them. When B is high enough (> 7T), all 
electrons in the emitter are in the lowest Landau level 
(n e = 0, m e < 0). These electrons can tunnel only into 
the QD states with the same non-positive angular mo- 
mentum. As a result when B > 7T tunnelling through 
the QD state rid = 0, rrid = 1 is forbidden. As a general 
selection rule, tunnelling through a state with angular 
momentum 777 is forbidden whenever there are only m 
occupied Landau subbands in the contacts. 

The figure [| shows that our theory is in good agree- 
ment with the experimental findings for the two first 
conductances peaks. However the complex theoretical 
behaviour of the third resonance is not observed. The 
theory predicts that the amplitude is given by the state 
7i = 0,777 = +1 below 6T and by 77 — 0, m = —2 at 
higher magnetic fields. Thus we expect to observe a de- 
crease of the current as B increases up to 6T, followed by 
a substantial increase and eventually a decrease at higher 
fields. Experimentally we observe only a decrease of the 
amplitude of this resonance. Nevertheless the coupling of 
the third and the fourth states (A and ▼) is clearly seen 
at 6T on the Darwin- Fock spectra (fig. IJd). Notice also 
that the right side slope of the third conductance peak 
(indicated by a bold line in fig.^a) is minimum at 6T. So 
we attribute the experimental decrease of the third reso- 
nance amplitude to a poor overlap with the first Landau 
level in the emitter, which decreases the contribution to 
the current of the n = 0,m = —2 state. Another expla- 
nation could be releated to the non parabolicity of the 
real QD. 

In conclusion, we have presented an analysis of the 
magneto-transport through QDs. We have shown that 
the current amplitude of the resonances is directly related 
to the angular momentum of the QD states contributing 
to the transport. The magnetic field strongly suppresses 
tunnelling through states with angular momentum paral- 
lel to it and in some cases tunnelling through such states 
is altogether forbidden. This analysis is substantiated by 
our experimental data collected on 3D-0D-3D structures 
and the agreement between the experimental findings and 
the theoretical results is very good. The analysis can 
be readily extended to 2D-0D-2D systems such as those 
where an accumulation layer is formed in the contacts || 
and, with more substantial modifications, to study QDs 
where the single particle picture no longer holds. 

We wish to express our gratitude to C. Mayeux and 
D. Arquey for their technical help. 
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FIG. 2. (a) I-V characteristics of the 35nm-radius at 
T — 35mK for different values of B. (b) Evolution of the 
I-V curve as B is increased from to 6T in steps of 0.2T. All 
the curves have been offset for clarity. 
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FIG. 3. (a) dl/dV spectra for several values of B from 0T 
to 15T (b) Voltage position of the dl/dV peaks as a function 
of B. The dotted lines represent the single-electron energy 
band diagram with a — 0.5 and hujo — 15meV. 
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FIG. 1. (a) Schema of the device, (b) Schematic energy 
band diagram of the device under an applied bias. A magnetic 
field is applied along the current direction, (c) Shift of the 
threshold voltage with LED illumination (1) and warm-up 
(2). After each LED pulse, the whole spectrum permanently 
shifts towards lower voltage bias. 
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FIG. 4. (a) Experimental current amplitude for the first 
three resonances of the 35nm-radius quantum dot device. The 
second and third resonances show, respectively, a maximum 
and a minimum at a magnetic field of about 6T. (b) Calcu- 
lated current amplitude as a function of B for different quan- 
tum dot states. Parts of these curves which cannot correspond 
to the experimental third resonance are dotted. Inset: Cur- 
rent contribution of the first Landau level in the emitter for 
n e = rid = and m — 0, — 1, —2 (from top to bottom). 
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